Effects of the electron-phonon coupling near and within the insulating Mott phase 
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The role of the electron-phonon interaction in the Holstein-Hubbard model is investigated in the 
metallic phase close to the Mott transition and in the insulating Mott phase. The model is studied 
by means of a variational slave boson technique. At half-filling, mean-field static quantities are in 
good agreement with the results obtained by numerical techniques. By taking into account gaussian 
fiuctuations, an analytic expression of the spectral density is derived in the Mott insulating phase 
showing that an increase of the electron-phonon coupling leads to a sensitive reduction of the Mott 
gap through a reduced effective repulsion. The relation of the results with recent experimental 
observations in strongly correlated systems is discussed. 
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I. INTRODUCTION 

The understanding of the properties of Mott insula- 
tors such as CaTiOs and ^203 represents a long-standing 
problem. As suggested by Mott^ the strong Coulomb re- 
pulsion among the electrons can cause a metal-insulator 
transition opening a gap in the density of states that 
is usually known as Mott gap. The correlation-driven 
metal-insulator transition is often investigated within the 
framework of the Hubbard lattice Hamiltonian which 
takes into account the electron interaction through the 
on-site repulsion term U. As developed by Hubbard;^ 
two subbands generically called lower and upper Hub- 
bard bands separated by the Mott gap of the order of 
the energy U arise in the excitation spectrum for large 
enough repulsion. 

The correlation effects are important not only on the 
insulating side but also on the metallic phase where a 
great deal of insight has been obtained by using the 
Gutzwiller wave function^ A progress toward the under- 
standing of the metal-insulator transition has been made 
by the reformulation of the problem in terms of slave 
bosons^ which reproduces the Gutzwiller approximation 
at the saddle-point level and allows to study the effects 
of gaussian fluctuations in the metallic^ and in the in- 
sulating phasei^ In this state the resulting single-particle 
spectral function consists of two broad incoherent con- 
tributions reminiscent of the lower and upper Hubbard 
bands. An overall description of the Mott transition in 
the Hubbard model can be obtained by means of the dy- 
namical mean field theory (DMFT) that is exact in the 
limit of infinite dimension.^ The study of the infinite di- 
mensional model has essentially confirmed the validity of 
the description obtained when slave boson fluctuations 
about the mean-field (MF) are considered. 

Although the Hubbard Hamiltonian captures the fun- 
damental properties of systems near the Mott transition, 
it does not take into account the role of the lattice degrees 
of freedom. Actually, the presence of strong electron- 
phonon (el — ph) interactions has been pointed out in 
several systems, such as cupratesjSiS manganitesjiS and 



^203-^^. In addition to the on-site repulsion, the most 
natural model incorporates a short-range interaction of 
the electrons with local phonon modes of constant fre- 
quency wqJ*^ This model is described by the Holstein- 
Hubbard Hamiltonian H 



t'^alai-\- guJo'^ni(^ai + al^ . (1) 
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Here t is the electron transfer integral between nearest 
neighbor sites < i,j > of a d-dimensional simple cubic 
lattice, (cio-) creates (destroys) an electron with spin 
a at the site i and Ui — riio- — ^l^Cia is the local 
density operator. In Eq.(^) a| (a^) is the creation (anni- 
hilation) phonon operator at the site i, and the parame- 
ter g represents the coupling constant between electrons 
and local displacements. The dimensionless parameter 
A = g^LUo/W, with W = 2dt bare electron half band- 
width, is typically used to measure the strength of the 
el — ph interaction in the adiabatic regime characterized 
by small values of the ratio 7 = coo/t. Within this regime 
the Coulomb repulsion is found to dominate the prop- 
erties of the metallic phase also with a sizable el — ph 
coupling. Furthermore, there is a very little softening 
of the phonon frequency on the approach to the Mott 
transition since the Hubbard term U suppresses charge 
fluctuationsii^ Actually within the Mott phase the spec- 
tral density shows phonon side bands whose peaks are 
separated by multiples of the bare frequency ujq^ Fi- 
nally, close to the Mott transition at finite density, an 
intermediate el — ph coupling leads to the phase separa- 
tion between a metal and an insulator. -'^^ 

In this paper the study of the Holstein-Hubbard model 
focuses on the role of the el — ph interaction in mod- 
ifying the physical properties of the electrons close to 
the metal-insulator transition and in the insulating Mott 
phase. The starting point of the approach is the gener- 
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alized Lang-Firsov transformation-'^^ U = e , with 

^ = f Xl[< > +fiini- < >)](ai - a|), (2) 

where the parameters ji take into account the polaronic 
local density fluctuations which couple to the lattice dis- 
tortions. The next step is the functional-integral repre- 
sentation in terms of the four-slave bosons e,, pio-, and di 
obtained by imposing the equivalence with the original 
model through the Lagrange multipliers a|^^ and A-^^'' 
First we consider the MF solution at half-filling, then 
slave boson gaussian fluctuations on the top of the MF 
finding that the resulting Mott gap is sensitively reduced 
with increasing A since it is determined by the el — el 
repulsion renormalized by the effects of the el — ph cou- 
pling. Finally the relation of the results with recent ex- 
perimental observations in strongly correlated systems is 
discussed. 

The MF solution in the paramagnetic homogeneous 
phase is obtained by replacing the Bose fields with their 
mean values ((ci) = eg, (di) = do, {pi„) = po) and by 
assuming /, = /, A^^ = a[,^\ and x[f = A[,^^ii The 
resulting mean-field Hamiltonian is characteristic of free 
fermion quasi-particles and phonons and the mean double 
occupation is controlled by Uefj = U + 2g'^ujQ{p ~ 2/). 

At half-filling the MF energy per site u is given by a 
simple functional depending on / and 

u = ge--^'s'e-HC/e//dg-g2cjo, (3) 

with q — 8(io ^ 16^0 and e mean bare kinetic energy. 
The Mott phase is the insulating state for large C/e/ / > 
characterized by g = 0, do = 0, and the energy per site 
equal to — g^wo (characteristic energy of the limit U jt = 
cx)). When C/e// becomes negative, a on-site bipolaron 
solution sets in corresponding to do — 0.5, / — 1, and 
energy u = C//2 — 2g^ujo. 

In Fig. 1 the phase diagram at half-filling is reported 
for 7 = 0.2 in the three-dimensional case. By analyzing 
the behavior of do , it is found that the transition to MI is 
always found to be second order, that to the bipolaronic 
insulator (BI) is first order, finally that between MI and 
BI is again discontinuous in excellent agreement with the 
results derived by the DMFT solution.^'^ In particular we 
notice that, with increasing A, the line separating the M 
and MI phases shifts to values of U larger than Uc, the 
critical value in absence of the el — ph coupling. Actually 
the Mott transition is influenced by the el — ph interac- 
tion since it is Ueff that governs the transition and it 
becomes smaller with increasing A. Therefore, the con- 
dition Ueff — Uc, characteristic of a transition driven by 
the el — el interaction, implies that the transition occurs 
for larger values of the bare U. Within the MF approach 
the interplay between el — el and el — ph interactions in 
affecting the Mott phase is essentially linked to the value 
of the parameter / that, at fixed adiabaticity ratio and 
A, is weakly decreasing with increasing U and it is of the 



order of 7/4d in the adiabatic regime near the Mott tran- 
sition. The transition line between Af and MI phase is 
given by A ~ {U-Uc)/{2fW) ~ {U~Uc)ht. Therefore, 
as shown in the inset of Fig. 1, the dependence of the 
Mott transition line on A becomes also more pronounced 
with increasing the adiabaticity ratio 7. In the atomic 
limit (7 = 00) we recovers the exact solution without 
metallic phase with / = 1 and Uc = Q (dotted line in 
Fig. 1 corresponding to A = U/2W). 

While the Mott transition is driven by the growth of 
the spin susceptibility, the transition from M to BI is 
characterized by the enhancement of the charge fluctu- 
ations inducing a decrease of the effective repulsion. In 
Fig. 2 we report the MF results for the spectral weight 
Z at the Fermi energy (equal to m/m*) and the local 
magnetic moment M =< {M^Y > stressing the effects 
of the el — ph interaction. Within the MF approach we 
simply have Z = qe~^ ^ and M = 1 — 2do. Far from the 
Mott transition {U /Uc small) Z decreases with increasing 
the el — ph coupling as expected for any localizing inter- 
action. However, near the Mott phase, the effects due 
to the reduction of Ueff become more marked and are 
able to induce the enhancement of Z with increasing A 
that, again, is in good agreement with DMFT resultsii^ 
In Fig. 2 we also show that the M phase is reduced in 
comparison with its value at A = for any ratio U /Uc 
implying that the double occupation do increases by ap- 
proaching the BI phase. Therefore, while M increases 
as function of C/, it decreases as function of A. 

The MF solution can be readily generalized at densi- 
ties different from half-filling. Within the MF approach 
some of us have shown that the interplay of strong local 
el— el and el—ph interactions can push the system toward 
a phase separation between states characterized by dif- 
ferent lattice distortions >ii The phase coexistence occurs 
for intermediate values of the el — ph coupling and its 
relevance within the Hubbard-Holstein model has been 
confirmed also by DM FT works:-^- 

The task of including charge fluctuations described by 
the e and d fields is simplified in the Mott phase by the 
fact that at MF level the Bose fields e and d are not 
condensed (eo = do = 0) , while pq ~ 1/ \/2& At gaussian 
level the fluctuations of the e and d fields are actually de- 
coupled not only from those of the p fields but also from 
the phononic a fields. Clearly the matrix elements of the 
fluctuating fields e and d are renormalized by el — ph 
contributions because of Uef / . The procedure of calcula- 
tion is the following: we start at densities different from 
half-filling (n < 1), then we calculate the quantities in 
the limit n ^ \^ for values of the parameters where the 
system is in the Mott phase at half-fillingi& 

The inclusion of gaussian fluctuations allows to cal- 
culate the one-particle spectral function A{lij) of the in- 
sulating Mott phase in the limit n ^ 1~ (at MF it is 
identically zero) . The incoherent contribution arises from 
the complicated motion of a quasi-particle surrounded by 
the cloud of charge and lattice excitations that it leaves 
behind. In the limit n ^ 1~ it is possible to derive an 
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analytic expression of the spectral function A{lu) yielding 
at zero temperature 

A{ij) = e-9'/'i(a,)+ (4) 

Tl=l 

^-9^f \ > 9{lu - nLUo)A{Lu - ncog), (5) 

•■^ — ' rt! 

Tl=l 

with 9{x) Heaviside function and the function A{uj) akin 
to that obtained for the Hubbard model with U = Ueff.^ 
The first term of Eq. ||SJ) is the product of two quanti- 
ties, with e~-^ ^ renormalization factor due to el — ph 
coupling. Through A{u!) this term is able to recon- 
struct lower and upper Hubbard bands. The second 
and the third term in Eq. ((SJ represent the contribu- 
tion due to the phonon replicas of hole- and particle- type, 
respectively,^^^ In Fig. 3(a) we report the spectral density 
N{u!) — A{uj)/2tt (solid line) together with the resulting 
first term of Eq. @ (dotted line) and the contribution 
from the phonon replicas (dashed line). We notice that 
this last term provides a non negligible spectral weight 
to the total spectral density at energies out of the gap 
of A{uj). Therefore, the Mott gap of the system is deter- 
mined by the gap of the function A(uj) and it is simply 
given by A = Ueff£,- As seen in Fig. 3(b), the reduc- 
tion of the gap with increasing A can be also of the order 
of bare half bandwidth Finally we stress that the 

calculated gap is traceable to the difference of the chem- 
ical potentials at n = 1 and n — 1^, in analogy with 
the results of the Hubbard modeli^ For n = 1 we have 
/Lt(l) = U/2 — 2g^(jJo, while for n = 1~ the chemical po- 
tential is /i(l^) — /i(l) — Ueff(,/2. Since the system has 
particle-hole symmetry, the gap A can be obtained as 

A = 2[Ml)-/.(l-)]=[/,//e 

In Fig. 4 we report the difference between the gap at 
finite A and that at A = as a function of the adiabatic- 
ity ratio. Since the attraction between the electrons gets 
larger with increasing 7, the resulting Mott gap is more 
reduced. However in the strong el—ph coupling regime (A 
larger than 1) there is no dependence on the adiabaticity 
ratio since the particles localized by the strong correla- 
tion are small polarons. In fact in this regime the MF 
solution at finite density is minimized for / = 1 yield- 
ing the effective interaction Ueff — U — 2g^uJo for any 
finite value of adiabatic ratio. In the inset of Fig. 4 we 
show the variation of the effective interaction Ueff as a 
function of the adiabaticity ratio making the comparison 
with the behavior of the Mott gap. In the regime where 
A < 1, the quantity ^ affects the magnitude of the gap 
and its dependence on the adiabaticity ratio, while in the 
strong coupling regime ^ ~ 1 implying that the gap in 
units of W is just Ueff/W = U/W - 2A. 

In this work we have seen that in the metallic phase 
close to the Mott transition the spectral weight Z is en- 
hanced and the gap in the insulating Mott phase is re- 
duced as A is raised. These behaviors can be related 



to some recent experimental and theoretical studies in 
V2O3 and Cr-doped ^203"^^°'^^ where the Mott gap is 
unexpectedly small and in the metallic phase the quasi- 
particle peak in the photoemission spectrum has a signif- 
icantly large weight in comparison with that theoretically 
predicted. We suggest that the inclusion of the el — ph 
interaction could be able to partially fill the discrepancies 
between the experimental observations and the theoreti- 
cal studies. Clearly the single orbital model is not suffi- 
cient to fully explain the electronic and magnetic struc- 
ture of such systemS)2S so a proper multi-orbital theory 
has to be considered in order to obtain a better agree- 
ment of the theory with experiments^ The results due 
to the reduction of the effective repulsion caused by the 
el — ph coupling are valid when the lattice distortions 
are coupled to charge fluctuations like in the model of 
Eq.Q. However for general models the issue is delicate 
since there are interactions such as the Jahn- Teller cou- 
pling for which the phonon-mediated attraction could be 
even diminished by strong correlations^ 



Finally we can argue the modifications of the spectral 
properties at densities near the MI phase and at half- 
filling just under the edge of localization. Clearly there is 
a finite spectral weight at the Fermi energy where quasi- 
particle states begin to appear. Therefore, in addition 
to the incoherent contribution at high energy, the spec- 
tral density presents also the coherent term. Close to the 
metal-insulator transition the coherent term could not be 
strongly affected by the el — ph coupling. In fact for the 
combined effect of the strong correlation and el — ph cou- 
pling the quasi-particle band can be narrower than the 
phonon energy luq implying the impossibility of the single 
phonon scattering between the quasi-particles. There- 
fore near the Mott transition the el — ph coupling af- 
fects mainly the incoherent term of the spectral density. 
This result is in agreement with recent experiments on 
Bi2Sr2CaCu20s+s made using angle-resolved photoe- 
mission spectroscopy.^ In fact it has been found that the 
oxygen isotope substitution mainly influences the broad 
high energy humps. 



In conclusion, we have discussed the role of the el — ph 
interaction in modifying the physical properties of the 
electrons in the metallic phase close to the Mott transi- 
tion and in the MI phase. The approach to study the 
Holstein-Hubbard model has been based on a variational 
slave boson technique that provides results in agreement 
with DM FT. An analytic expression of the spectral den- 
sity is derived in the Mott phase showing that due to the 
reduced effective repulsion the Mott gap decreases as the 
el — ph coupling constant increases. In this paper we 
have mainly discussed the phases without long-range or- 
der. The study of broken symmetry phases is possible 
within the slave-boson formalism^ and it is left for future 
work. 
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Figure captions 

Fl The phase diagram U/W versus A at half-filhng 
for 7 = 0.2 in the three-dimensional case. The 
transition hnes separate the metaUic state M from 
the Mott insulator MI and the bipolaronic in- 
sulator BI. The dotted line is the locus where 
U — 2g'^oJo = 0. In the inset, the transition line 
between the metallic and Mott insulating phase is 
shown for different adiabaticity ratios 7. 

F2 The difference between the spectral weight Z at the 
Fermi energy (a) and the local magnetic moment 
M at finite A (b) and their respective values at 
A = for some values of the ratio U/Uc in the 
three-dimensional simple cubic lattice. 

F3 (a) The renormalized density of states N (solid line) 



2 j-2 

together with the dominant contribution ^ 
A{uj) (dotted line) and the term due to the phonon 
replicas (dashed line) as function of the frequency 

UJ. 

(b) The renormalized density of states TV for differ- 
ent values of the el — ph coupling constant A as a 
function of the frequency w. 

F4 The difference between the gap at A = and that at 
finite A in units of the bare half bandwidth W as 
function of the adiabaticity ratio for different values 
of the el — ph coupling. In the inset the effective 
repulsion Uef f and the Mott gap A as a function 
of the adiabaticity ratio 7. 
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